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Abstract 

We consider n non-intersecting Brownian motion paths with p prescribed starting po- 
sitions at time t — and q prescribed ending positions at time t = 1. The positions of 
the paths at any intermediate time are a determinantal point process, which in the case 
p = 1 is equivalent to the eigenvalue distribution of a random matrix from the Gaussian 
unitary ensemble with external source. For general p and q, we show that if a temperature 
parameter is sufficiently small, then the distribution of the Brownian paths is characterized 
in the large n limit by a vector equilibrium problem with an interaction matrix that is 
based on a bipartite planar graph. Our proof is based on a steepest descent analysis of an 
associated (p + q) x (p + q) matrix valued Riemann-Hilbert problem whose solution is built 
out of multiple orthogonal polynomials. A new feature of the steepest descent analysis is a 
systematic opening of a large number of global lenses. 

Keywords: non-intersecting Brownian motions, Karlin-McGregor theorem, vector po- 
tential theory, graph theory, multiple orthogonal polynomials, Riemann-Hilbert problem, 
Deift-Zhou steepest descent analysis. 

1 Introduction 

This paper deals with non-intersecting one-dimensional Brownian motions with prescribed start- 
ing and ending positions. This model has already been discussed in various regimes. For the 
case of one starting point and one ending point it is known that the positions of the paths at 
any intermediate time have the same distribution (up to trivial scaling) as the eigenvalues of a 
Gaussian Unitary Ensemble from random matrix theory |18j . Moreover, as the number of paths 
tends to infinity and after appropriate scaling, the paths fill out an ellipse in the ix-plane, see 
Figure [TJ 

In the case of one starting point and two or more ending points the positions of the paths 
have the same distribution as the eigenvalues of a Gaussian Unitary Ensemble with external 
source. This model is described by multiple Hermite polynomials. As the number of paths tends 
to infinity, the paths fill out a more complicated region whose boundary has cusp points. The 
limiting distributions can be computed in terms of an algebraic curve known as Pastur's equation 
[21 03 [2H HH1 See Figure [2] for an illustration of the case of two ending points. 



*Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 200B, B-3001 Leuven, Belgium, 
email: {steven. delvaux, arno. kuijlaars}@wis. kuleuven.be. 

The first author is a Postdoctoral Fellow of the Fund for Scientific Research - Flanders (Belgium). 
The work of the second author is supported by FWO-Flanders project G. 0427. 09, by K.U. Leuven research grant 
OT/08/33, by the Belgian Interuniversity Attraction Pole P06/02, by the European Science Foundation Program 
MISGAM, and by grant MTM2008-06689-C02-01 of the Spanish Ministry of Science and Innovation. 



1 



-1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 1: Non-intersecting Brownian motions with one starting and one ending point. As the 
number of paths tends to infinity, the paths fill out an ellipse in the ix-plane. 
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Figure 2: Non- intersecting Brownian motions with one starting and two ending positions. 

The case of one ending point and two or more starting points is equivalent due to the time 
reversal symmetry in the model. 

Much less is known for the general case of p > 2 starting points and q > 2 ending points. 
For example, it is not known whether there exists an underlying random matrix model for this 
case. What is known is that the model is described by multiple Hermite polynomials of mixed 
type [5] which have a characterization in terms of a (p + q) x (p + q) matrix valued Riemann 
Hilbert problem. Calculations for the limiting distributions of paths in the large n limit were 
done for very specific cases with p = q = 2 in [TUJ \W\ based on the spectral curve (analogue of 
the Pastur's equation) that could be computed in these special cases, see also the related works 

Dug. 

It is the goal of this paper to study the case of general p > 2 and q > 2 with methods from 
potential theory, more precisely with vector equilibrium problems with external fields. 

Equilibrium problems with external fields were developed in approximation theory in the 
context of orthogonal polynomials, Padc approximation, and polynomial approximation with 
varying weights, see |30| I31j . They are also a powerful tool in the study of unitary random 
matrix ensembles [TTJ [T2] . Vector equilibrium problems were studied in the context of Hermitc- 
Pade approximation and the associated multiple orthogonal polynomials [2J [57J . 
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Figure 3: Non-intersecting Brownian motions with two starting and two ending positions. The 
starting and ending positions are sufficiently far apart so that around time t — 1/2, there are 
three groups of paths in the large n limit. 

For p = q = 2 we are considering a situation such as the one shown in Figure [31 where a 
certain fraction of the paths starts in each of the two starting points, and ends at each of the two 
ending points. As the number of paths increases we see the following situation. For small time 
the paths are in two separate groups that emanate from the two starting positions. At a certain 
time one of the groups splits into two, leading to a situation of three separate groups of paths. 
Then at a later time two of the groups come together and we end up with two groups that end 
at the two ending points. 

Our results will deal with the situation at times where there are three groups of paths, or for 
general p and q, where there are the maximal number (namely p + q — 1) of groups of paths. 

There is an alternative possible scenario in which the two groups of paths first merge into 
one group and later split again into two groups of paths. This will happen if the starting and 
ending positions are sufficiently close to each other. The first scenario happens if the starting and 
ending positions are relatively far away from each other. Below we will actually distinguish the 
two scenarios in terms of a temperature parameter T so that for small T we have the situation 
with the three groups of paths. 

2 Statement of results 
2.1 Assumptions 

Let p > 2 and q > 2. We fix p starting points oi, 

Qi > a 2 > 

and q ending points b± , . . . , b q with 

h > b 2 > ■ ■ ■ > b q . (2.2) 

For a given (large) n we consider n non-intersecting Brownian motion paths and we assume 
that rik of the paths start at Ofc and that mi of the paths end at 6/ for k = 1, . . . ,p and I = 1, . . . , q. 
Thus 

p q 
y^n fc = y^m; = n. 

fc=i i=i 



. , a p which we assume to be ordered as 
> a p , (2.1) 
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Since the paths are non-intersecting, the numbers n k and mi also determine for each k = 1, . . . ,p 
and I = 1, . . . , q, the number n kt i of paths that start at a k and end at b\. We call the fractions 

4? = ^ (2-3) 

the finite n transition numbers. Note that 

4?>o, EE4" ) = i- (2.4) 

As n — > oo, we assume that the finite n transition numbers have limits 

t k .i = lim 4™ } (2.5) 

which are the limiting transition numbers. It is convenient to arrange the (finite n and limiting) 
transition numbers into p x q matrices 



, . ' (* fe .')fc=i...., P ,i=i. 



To avoid degenerate cases, we assume that each row and column of the matrix (tfe,i)fe = i p ; = i q 
has at least one non-zero entry. 

The assumption that the paths are non-intersecting puts a number of constraints on the 
numbers nk,i and on the limiting transition numbers tk t i- Indeed, not all au can be connected to 
all bi and certain transition numbers must be zero. The constraints on the transition numbers 
are easy to visualize in terms of a weighted bipartite graph 

G=(V,E,t), (2.6) 

with vertices 

V = {oi, . . . , a p } U . . . , bq\, (disjoint union), 

edges 

E = {(a k M) G V x V | i fci , > 0}, 

and a weight function 

t (0,1] : (afcjfcO^tfc,!- (2-7) 

Example 1. TTie graph G — (V,E,t) associated with Figured is shown in Figure^ The graph 
has four vertices and three edges, each of them with weight 1/3. 




Figure 4: The graph associated with Figure [31 
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Figure 5: Non-intersecting Brownian motions with two starting points and four ending points, 
as discussed in Example [3 

Example 2. For a more complicated example we consider a situation with p — 2 starting points 
and q = 4 ending points as in Figure [21 
The matrix of transition numbers is 



The constraints on the transition numbers are contained in the following obvious result that 
we state without proof. 

Proposition 2.1. The graph G has the following properties: 

(a) G has at most p+q—l edges. For each i = 1, . . . ,p+q—l, there is at most one non-vanishing 
transition number tk,i with k + I — 1 = %■ 

(b) G is a connected graph if and only if the number of edges is equal to p + q — 1 . 

(c) G has no cycles (and so G is a tree if G is connected). 
In [TOl [H] the special case of transition numbers 




(2.8) 




Figure 6: The graph associated with the transition numbers f[2.8}> - 
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was considered. This example leads to a non-connected graph. 

In the case of a connected graph the structure of the matrix of transition numbers (tk,i) is 
easy to describe. 

Proposition 2.2. Suppose that the graph G is connected. Then the non-zero entries of the 
matrix \tk,i) k i are situated on a right-down path starting at the top left entry (1, 1) and ending 
at the bottom right entry (p, q) . The steps in the path are either by one unit to the right ( a right 
step) or one unit down (a down step). 

In Examples [T] and [5] we have 

fx 0\ fx x 0\ 

\x x) ' ^0 x x x) ' 

respectively. 

In this paper we consider only connected graphs, and we will make the following assumption. 

Assumption 2.3. We assume that the graph G is connected. That is, we assume that G has 
p + q — 1 edges, and for every i = X, . . . ,p + q — 1 there is exactly one non-vanishing transition 
number tk.i with k + I — 1 = i, and we define 

k(i) = k, l{i) = I, if tk,i > 0, and k + I — 1 = i. 

It follows from the assumption and from (|2.5|) that also for large enough n, the finite n 

(n) 

transition numbers have the same non-zero pattern. Thus t k ; > if and only if tk,i > 0. 

In what follows we follow the convention that i labels the edges E of the graph, and so we 
identify i with the edge 

( a k(i) i 

of the graph. 



whose overall variance 



2.2 Non-intersecting Brownian motions 

We consider Brownian motions having transition probability density 

P(t, x, y) = -^L=e-^-^ (2.9) 

a 2 = - (2.10) 

n 

is proportional to 1/n where n is the number of Brownian paths. We interpret the proportionality 
constant T > as a temperature variable. 

In this paper we consider small temperature T. We show that for small T the paths at time 
t have a limiting mean distribution that is characterized by a vector equilibrium problem. In the 
first theorem we state the existence of a limiting mean distribution. 

Theorem 2.4. Consider n independent Brownian motions with transition probability 
(|2.10p conditioned so that 

• the paths are non-intersecting in the time interval (0, 1), 

• nk.i of the paths start at and end at bi, for each k = 1, . . . ,p and I = 1, . . . , q. 



G 



Assume that as n — > oo, the finite n transition numbers nk^i/n converge to tk,i, and that the 
corresponding graph G is connected. 

Let t £ (0, 1). Then there exists a T* = T*(t) > so that for all T £ (0, T*) the limiting 
mean distribution of the positions of the paths at time t exists, and is supported on the union of 
p + q — 1 disjoint intervals US 7-1 [ a ii Pi\ with a density pi on the ith interval. 

The vector of measures (fix, . . . where dfj,i(x) = pi(x)dx for i = 1, . . . ,p + q — 1 is 

the minimizer of a vector equilibrium problem that will be described in the next subsection. 



The proof of Theorem 12.41 will be based on a Deift-Zhou steepest descent analysis of the 
Riemann-Hilbert problem in Section 12.51 The details of the steepest descent analysis will be 
described in Sections [M3 and the proof of Theorem 12.41 will then be given in Section [51 

Remark 2.5. The special case where p = 1, q = 2 and mi = m,2 = n/2 was studied in (SKUEj. 
In [5] the energy functional (|2.15j) is given for this special case, but it was not used in the further 
analysis. Instead the results in [5] were stated and proved in terms of an algebraic curve (Pastur's 
equation). 

Remark 2.6. As already mentioned, for the case p = 1 and q > 2 the non-intersecting Brownian 
motions are distributed like the eigenvalues of the Gaussian unitary ensemble with external 
source. 

There are more general random matrix ensembles with external source, which however do 
not have an equivalent interpretation in terms of non-intersecting paths. The case of a quartic 
potential was studied in 25j. An analogue of Theorem 12.41 is valid in this case as well, provided 
that a suitable analogue of the temperature T is sufficiently small. 

Note that Theorem 12 .41 gives only a result for T sufficiently small, but it does not specify how 
small T should be. We expect that the theorem remains valid for all temperatures T that are 
such that at time t we have the maximal number of groups of paths. At a critical temperature 
T C Ht(t) we expect that two (or maybe more) neighboring intervals [a^/Jj] and [atj-fi, Ai+i] merge 
and a Pearcey phase transition occurs. However we were unable to prove this. 

2.3 Equilibrium problem 

The logarithmic energy of a measure on R is defined as usual by 

/(/i) = //log dn(x) dfi(y). (2.11) 

J J \x-y\ 

The mutual energy of two measures fi, v is defined by 

/(/x, v) = ff log dfx(x) dv{y). (2.12) 

J J \x-y\ 

We write for i = 1, . . . , p + q — 1, 

Xi{t) = (l-t)a Hi) +tb m , 0<t<l (2.13) 
and given t 6 (0, 1) we define the quadratic functions 

Vi(x)= — i— (s-sjft)) 2 , xeR, (2.14) 

for i = 1, . . . ,p+q— 1. These functions will play the role of external fields in the vector equilibrium 
problem that is relevant for our problem. 
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Definition 2.7. Fix t € (0, 1) and let T > 0. Consider the energy functional 



p+g-i * p+g-i „ 

E(fii,...,fx p+g -i) = a i,j I (Mi,fJ'j)+j; Vi(x)dfXi(x), (2.15) 

where is defined in (|2.14p . and the interaction matrix A — (ffljj) has entries 

1 if i = j, 

, ^ i if i ? j and fc(t) = fc(j) or Z(t) = l(j), (2.16) 
otherwise. 

Then the vector equilibrium problem consists in minimizing the energy functional (|2.15|) over all 
vectors of positive measures (nit . . . , Hp+q-i) supported on the real line for which 

dUi = *k(t),H(t)> for i = 1, . . . ,p + g - 1. (2.17) 

One may understand the energy minimization problem in Definition 12.71 as follows. To each 
of the edges of the graph G = (V, E, t) we associate a measure pn, i = 1, . . . ,p + q — 1, of total 
mass equal to the weight twjym of that edge. This measure represents a distribution of charged 
particles on the real line that repel each other due to the diagonal term 

in (|2.15p . For the particles of different measures Hi, Hj> * h there are two possibilities. The 
first case is when the entry of (|2.16p equals 1/2. This happens if the edges corresponding 
to i and j are adjacent in the graph G. Then there is repulsion between the measures Hi an( i 
Hj but with a strength that is only half as strong as the repulsion for each individual measure. 
The second case is when the (i,j) entry of (|2.16p equals zero. In that case there is no direct 
interaction between the measures Hi an( i Hj- 

The last term of (|2.15p is a sum of external field terms due to the action of the external 
field ikVi{x) on the measure Hi- The energy minimizer (hi, ■ ■ ■ ,Hp+q-i) m Definition 12.71 then 
corresponds to the equilibrium distribution of charged particles under the energy functional 

Proposition 2.8. The interaction matrix A is positive definite. 
Proof. It is easy to check that the interaction matrix A is equal to 

A=^B T B (2.18) 

where B is the incidence matrix of the graph G. That is, we choose a numbering k = 1, . . . ,p+ q 
of the vertices, and then we have 

B = (bk,i)k=l,...,p+q,i—l,...,p+q-l 

where bk,i — 1 if vertex k is incident to edge i, and otherwise. 

From (|2.18p we get that A is positive semi-definite, and for any column vector x of length 
p + q — 1 we have 

x T ^lx= i||Sx|| 2 > 0. (2.19) 
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Now assume that £?x = 0. Consider a leaf of G, i.e., a vertex which is incident to exactly one 
edge. Then B has exactly one zero in the row corresponding to this vertex, and from 6x = it 
follows that the component of x corresponding to the edge that is incident to the leaf vanishes. 
Since G is a tree (see Proposition I2.1f c)) we can then gradually undress G by peeling off leafs 
one by one and we conclude in this way that all components of x are equal to 0. Thus x = if 
Bx = 0, which implies in view of (|2.19[) that A is positive definite. □ 

Corollary 2.9. The vector equilibrium problem of Definition \2.l\ has a unique solution (pi, . . . , p, p + q -i) 
and each measure pi is compactly supported. 

Proof. The interaction matrix is positive definite by Proposition 12.81 The external fields Vi in 
the energy functional (|2.15|) have enough increase at ±oo so that standard arguments of potential 
theory as in [111 1271 130j can be used to establish the existence and uniqueness of the minimizer 
as well as the fact that each measure fii is supported on a compact set. □ 

Our next theorem describes the structure of the solution of the vector equilibrium problem 
for small T. 

Theorem 2.10. Fixt € (0, 1), and let (tfe,;) be a matrix of transition numbers. Then there exists 
T* > (the same T* that makes Theorem \2.4\ work) so that for every T 6 (0,T*) the following 
holds. 

(a) Each pi is supported on an interval 

supp(^) = [ai,[3i], i = l,...,p + q-l. 
The intervals [an , Pi] are pairwise disjoint and satisfy 

Pi+i <ati, i = l,...,p + q-2. (2.20) 

(b) The measure pn has a density pi with respect to Lebesgue measure which is real analytic 
and positive in the open interval (aii,/3i) and vanishes like a square root at the endpoints 
of [cii,f3i], i.e., there exist non-zero constants p^p and pf^ such that 

Pl {x)=pf ) ^/x~~a~+0{{x -a.fi 2 ) as x J. on, (2.21) 
Pi {x) = pi 2) VA -x + 0(([3i - xf' 2 ) asx^fSi. (2.22) 

The proof of Theorem 12.101 will be given in Section 13.11 
2.4 Special cases 

In Subsections 12.4. 1T|2. 4. 31 our main Theorem 12.41 will be illustrated for some special cases. 
2.4.1 The case p = 1: Angelesco-type interaction 

For the case p = 1 of one starting point, and an arbitrary number q of ending points, the graph 
G has a single vertex a\ on the left which is connected to each of the vertices b\ , . . . , b q on the 
right. For example, if q = 3 the graph has the form shown in Figure [71 
The energy functional (|2.15|) - (|2.16p is then equal to 

E(m,..., p.p+q-1 ) := I (^i) + 2^2 I (l J ' i '^ + T 51 j Vi{x)dpi{x). (2.23) 

i=l ijtj i=l ^ 
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Figure 7: A graph G with p = 1 starting points and q — 3 ending points. 



The functional (|2.23|) is exactly the one familiar from the theory of Angelesco systems [21] . All 
off-diagonal entries in the interaction matrix A in (|2.16p are equal to 1/2. For the example in 
Figure [7] the interaction matrix is 



1 1/2 1/2^ 
.1=1 1/2 1 1/2 
,1/2 1/2 1 



(2.24) 



2.4.2 The 'zigzag' case: nearest neighbor interaction 

Next we consider the case where p = q and the corresponding lattice path follows a zigzag line. 
The graph G is then just a chain of vertices: see Figure [HI 




a 2 b 2 a 3 



(b) 



Figure 8: For p = q = 3, Figure 8(a) shows a graph G which has zigzag form. Figure [8 (b) | shows 
the same graph written as a chain. 



The energy functional (|2.15[) now takes the form 

p+q— 1 p+q—2 



p+q-1 



i):= ^ J( M i)+ J! ^(W.W+iJ + y ^ Vi(x)dfH{x). (2.25) 
i=l t=l i=l 



We see that the interaction matrix of (|2.25p is tridiagonal with diagonal entries equal to 1, 
and entries on the first sub- and superdiagonal equal to 1/2. For the example in Figure [HI the 
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interaction matrix equals 



A = 



( 1 

1/2 




V o 



1/2 
1 

1/2 







1/2 
1 

1/2 







1/2 

1 
1/2 






1/2 
1 



(2.26) 



Note that the tridiagonal structure of the interaction matrix means that there is only nearest 
neighbor interaction. The neighboring measures repel each other, since all signs in the interaction 
matrix are positive. 

In a Nikishin system the interaction matrix is also tridiagonal, but the off-diagonal entries 
are —1/2 instead of 1/2, see |21j . 



2.4.3 The general case: block nearest neighbor interaction with Angelesco-type 
blocks 

Finally we consider the graph G in Figure El In this case the interaction matrix equals 



A 



( 1 

1/2 



V o 



1/2 
1 

1/2 






1/2 
1 

1/2 
1/2 





1/2 
1 

1/2 



\ 


1/2 
1/2 

1 / 



(2.27) 



Note that this interaction matrix is a mixture of the nearest neighbor and Angelesco interaction 
matrices. More precisely one could say that (|2.27|) has 'block' nearest neighbor interaction, 
where each of the blocks in turn has an Angelesco-type interaction, and with subsequent blocks 
intersecting in exactly one entry. For the matrix in (|2.27p these building blocks are 



1 

1/2 



1/2 
1 



1 

1/2 



1/2 
1 




2.4.4 Comparison with [21j 

A graph-based vector equilibrium problem was also considered in the work of Gonchar, Rakhmanov, 
and Sorokin |21j . The rule for building the interaction matrix from a graph G is similar to the 
one in this paper. The vector equilibrium problem is also labeled by the edges of a graph, which 
in [21], however, is a directed rooted tree. The off-diagonal entries of the interaction matrix are 
non-zero if the corresponding two edges have a common vertex. The entry is 1/2 if the two edges 
have a common initial vertex, and —1/2 if the final vertex of one edge agrees with the initial 
vector of the other. 

In our case, we only have 1/2 since the latter situation cannot happen. 



2.5 About the proof of Theorem 12.41 

The proof of Theorem 12.41 uses the connection of non-intersecting Brownian motions with pre- 
scribed starting and ending points with a determinantal point process and an associated Riemann- 
Hilbert problem. We recall this connection. 
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2.5.1 Determinantal point process 

The positions at time t € (0, 1) of n non-intersecting Brownian motions, starting at distinct dj, 
j = 1, ... , n., and ending at distinct positions bj, j = 1, . . . , n, have the joint probability density 
function 

- det (P(t, ai,Xj))" j=1 • det (P(l - i, Xi, &i))"j=i , 

with the transition probability density P defined in (|2.9|) and with Z a normalization constant. 
This is a consequence of a theorem of Karlin and McGregor [23] . (For applications of the discrete 
version of the Karlin-McGregor theorem see e.g. [22] .) In the confluent limit where n k of the 
starting positions come together at a k , k = 1, . . . ,p, and mi of the ending positions come together 
at bi for I = 1, . . . , q, the joint p.d.f. for the positions of the paths at time t can be written as 

V( Xl ,...,x n ) = ^det(/ i (x i ))" J=1 -det( 5i (x J -))" j=1 , (2.28) 
for certain functions fi, gi, that are built out of the p + q functions 

w 1>k (x) = e-^ (x - ak) \ k = l ) ...,p, (2.29) 
w 2 ,i(x) = e~^^ (x ~ bl) \ l = l,...,q, (2.30) 

see e.g. [5]- 

The p.d.f. (|2.28p defines a determinantal point process (in fact a biorthogonal ensemble, see 
[7]) with a correlation kernel K(x,y) that is such that 

P(xi,...,x„) = det (if(a:i,Xj))" i=1 (2.31) 

and for each m = 1, . . . , n, 

V(xi, ...,x n ) dx m+ i ■ ■ -dx n = ^ J"^' det {K(x u x ))™ j=1 . (2.32) 
In particular for m = 1, we have that 

—K(x, x) 
n 

is the mean density of paths. 

2.5.2 Riemann-Hilbert problem 

The kernel K can be described in terms of the following Riemann-Hilbert problem (RH problem) 
introduced in [9]. 

RH problem 2.11. The RH problem consists in finding a matrix-valued function Y : C \ R — > 

C (p+«)x(p+ ? ) such that 

(1) Y is analytic in C \ M; 

(2) For x € K, t/ie limiting values 

Y+(x) = lim F(z), F_(x) = lim F(z) 

z—*x, Im 2:>0 ^ — KE, Im ^<0 
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exist and satisfy 



Y + {x)=Y_{x) 



I P W(x) 
I„ 



(2.33) 



where Ik denotes the identity matrix of size k, and where W(x) denotes the rank-one matrix 
( outer product of two vectors ) 



W(x) 



(w 2}l (x) ••■ w 2tq [x)) (2.34) 

\lVl.pix); 

with w\ t k{x), k — l,...,p, and w 2 ,i(x), I = 1, . . . , q given by (|2.29p and (|2.30p . 
(3) As z — > oo, we have that 

Y(z) = (I p+q + 0(1/ z)) diag (z" 1 , . . . , z n ",z~ m \ . . . , z- ro «) . (2.35) 

The RH problem has a unique solution that can be described in terms of multiple Hcrmitc 
polynomials. This is shown in [9], generalizing results in [T9l [32j . According to [9j the correlation 
kernel is expressed in terms of the solution to the RH problem as 



K(x,y) 



2iri(x — y) 



(0 ••• w 2>1 (y) ... w %q {y)) Y~ 1 (y)Y + (x) 



fwi tl (x)\ 



w l,p( x ) 




V J 



(2.36) 



It is also worth noticing that if Yi^z) denotes the top leftmost p by p block of the RH matrix 
Y(z), then the determinant of ij. i(z) equals the average characteristic polynomial 



detY M (.z) = E 



Y[(z Xj 

3=1 



where the expectation E is taken according to the joint probability density ()2.28|) . This was 
shown in [4] for the case p = 1 and in [15j for the case of general p and q. 

2.5.3 Asymptotic analysis 

We analyze the RH problem 12.111 in the limit described in Theorem 12.41 That is, we take 
a 1 = T/n and we let —>■ oo, mi —>■ oo, so that 



Tiki 



tk,i 



asn^ oo. 



If, for each n, we denote the correlation kernel (|2.36jl by K n , then the limiting density of 
paths is 

p(x)= lim -KJx,x), i£i (2.37) 
n— *oo n 

The proof of Theorem 12 . 41 will be based on a steepest descent analysis of the RH problem l2.11l 
As a byproduct of this analysis we can also show that the local scaling limits of the limiting 
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distribution of the Brownian motions are those familiar from random matrix theory, i.e., they 
are described in terms of the sine kernel in the bulk and the Airy kernel at the edge. We will 
not discuss this any further and refer to the papers El EH HE] for a more detailed analysis in 
a similar context. 

2.6 Outline of the rest of the paper 

The rest of this paper is organized as follows. In Section [3] we establish general properties of the 
equilibrium problem, in particular leading to the proof of Theorem 12.101 In Section |4] we define 
the ^-functions, A-functions and the associated Riemann surface. These functions are used in 
Section[5]to normalize the RH problem at infinity. In Section[5]we apply Gaussian elimination to 
the RH problem by opening global lenses, thereby making the RH problem locally of size 2x2. 
This construction is fully systematic and may have interest in its own right. The remaining steps 
of the RH steepest descent analysis are described in Section [7] Finally, in Section [8] we prove the 
main Theorem 12.41 

3 Properties of the equilibrium measures 

In this section we study the vector equilibrium problem of Definition 12.71 In particular we prove 
Theorem \TM 

3.1 Proof of Theorem ETTOl 

The main difficulty is in the proof of part (a) of Theorem 12.101 For this we use the following 
lemma. Recall that t £ (0, 1) is fixed. 

Lemma 3.1. For every e > and t > 0, there exists T £ > so that for every T € (0,T £ ) the 
following holds. 

If(jj,i, . . . ,^p-|_ g _i) is the minimizer for the energy functional (|2.15p under the normalization 
(|2~T7|) with 

min**(i),n(i) > r > 0, (3.1) 

then the support of /li is contained in [xi(t) — e, Xi(t) + e] for every i — 1, . . . ,p + q — 1. 
Proof. We are going to prove the following upper and lower bounds for the quantity 

E* := E(ni, . . . ,/ip+,_i) 
(which depends on the chosen normalization (|2.17l) ). 

(a) There exist constants C\ > and 62,5 such that for every S > and T > we have 

E*<C 2 j + C^. (3.2) 
The constant C± is independent of S and T, and C^a is independent of T. 

(b) There exist positive constants C3, and C4, such that for every T £ (0, 1] and every i = 
1, . . . ,p + q — 1, we have 

E* > -C 3 + ^(x- x t (t)) 2 , for x £ aupp(/ij). (3.3) 

In addition, all constants Ci, 62,5, C3 and C4 can be taken independently of the normalization 
(|2.17|) . and only C4 depends on r. 
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Proof of (a): We may assume that <5 > is sufficiently small so that 

mxa(xi(t) - x i+ i(t)) > 25. 

i 

Then the intervals [xi{t) — 5,Xi(t) + 8] are mutually disjoint. 

We fix transition numbers t^u^u^ > so that X}i^Wi),Z(») = ^' ^ n eacn °f the intervals 
[xi(t) — 5,Xi(t) + 6] we choose a measure Xi of finite energy with total mass 1 (for example, 
an appropriately rescaled and centered semi-circle law will do) and we put Vi = tku\ui\\i- We 
choose the measures Vi independent of T and so 

p+q-l 

a hJ I ( v '^ l 'j) ( 3 - 4 ) 
does not depend on T. The sum ()3.4[) depends continuously on the numbers ifc(j),/(i) and so 

p+q-l 

G-ifi — max ai t jI(vi,Pj) (3.5) 

exists, where the maximum is taken over all choices tk(i),i(i) > with ^ i tk(i),l(i) = 1- 
Since Vi (x) < 2 t(i-t) on t Xi (*) — ^' x » (*) + ^] j we a ^ so nave 

V^(a;)di/i(a;) < 



2i(l - t) 



and 

p+q-l 



^ / V t (x) dis^x) < d6 2 (3.6) 
i=l ^ 

with a constant C\ that is independent of T and 5, and also of the ife(jv;(j). Then part (a) follows 
from (|3.5p . (|3.6p , and the fact that 



p+a— l p+g-i „ 

E* < E(vi,...,v p+q -i) = Y OijI (v^ Uj) + — Y Vi(x)dvi(x). 



Proof of (b): For each i we take a measure Vi of total mass \vi\ = \\ni\\ and satisfying for 
some constant K > 1, 

Vi<Kfj,i, for i = 1, . . . ,p + q< — 1. (3-7) 
Then (1 — A)/i, + Az/; is a positive measure for every A S [— 1], and so 

E* < E ((1 - A)/ii + Ai*, . . . , (1 - A)/ip+,_i + Ai/p+,_i) (3.8) 

for every A S [— 1/iiT, 1], with equality for A = 0. Thus the A-derivative of the right-hand side of 
vanishes for A = 0, which leads to 



p+q-l p+q-l „ 

E* = a^I(^,^) + — I Vi(x)(diH(x)+dvi(x)). (3.9) 
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From the elementary inequality \x — y\ < \/x 2 + l\/y 2 + 1 it follows that 
I(jH,"j) = J J log j^-^dfii(x)dvj(y) 

>~ if ^g(x 2 + l)dni(x)d^(y) -- \og{y 2 + l)d^{x)du 3 {y) 



1 

> — 
- 2 



log(a; 2 + 1) {dUiix) + duj(x)) 



where for the last inequality we used the facts that log(a; 2 + 1) > and \\vj\\ < 1, < 1. 

Since all a,j > it then follows from (|3.9p that 



p+q— i /. , p+q— l 



E* > -- ^ a i>3 - / log(a; 2 + 1) (cZ^(x) + c^-(x)) + — 51/ (^0*0 + ^(a;)) 

i,i=l ^ i=l 

p+g-1 . / P+9-1 \ 

= 2T J W" 7 E ^log(^ 2 + l)j (rfMiW+^iCx)). (3.10) 

For i = 1, . . . ,£>+(? — 1 and T < 1, the functions Vi(a;) — T fog^ 2 + 1) are bounded 

from below and assume their minimum value in a fixed compact interval (independent of i and 
T < 1). This is clear from the explicit form (|2.14jl of the functions Vt, and from the fact that 
Y^tr 1 <Hj < H 1 for each i. Thus 

Vi{x) - T J2 a id fog^ 2 + !) ^ ~ C ° T " ~ C ° T ' ,T G M ' 

i=i 

for some constant Co > independent of T < 1 and i. 

Using this in (|3 . 1 Op we find that for each i = 1, . . . ,p + q — 1, 

E*>-C 3 + -^f Vi(x) dvi(x), (3.11) 

where C3 > is a constant independent of i and T < 1. The inequality (13 . 1 1 [) (with the same 
constant C3) holds for all measures with ||^|| = ifc(i).;(z) an d satisfying (|3.7p for some if. For 
every x in the support of m, we can approximate the point mass 5 X by such z^. It thus 
follows from (|3 . 1 1 1) that 

E* > -C 3 + -^=Vi{x) t k (i),i(i), for x € supp(/ij). 

If tk(i),i(i) > t as in (|3.ip . then we obtain (|3.3p with 

1 1 



C 4 



4 2t(l-i) 



Conclusion of the proof: Combining ()3.2|1 and (|3.3p we find that there exist positive con- 
stants C5 = C1/C4 and (76,5 = (C^s + C-$)jC± so that 

(a; — Xi (<)) 2 < C* 5 <5 2 + C 6 , 5 T, for x E supp(^), 
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for every i = 1, . . . ,p + q — 1, T < 1 and 5 > 0. 

Let e > be given. Choose first S > so that C$5 2 < \e 2 and then choose T £ £ (0, 1] so that 
C 6 ,sT £ < \e 2 . Then for every i = 1, . . . ,p + q - 1 and T < T £ , 

(x — Xi{t)) 2 < e 2 , for x € supp(/ii), 
and Lemma EUJ follows. □ 
Having Lemma 13. II the proof of Theorem 12. 101 is rather straightforward. 

Proof. (Proof of Theorem [2TT0|) 

(a) Let e > be such that the intervals [xi(t) — e, Xi(t) + e], i = 1, . . . , p + q — 1, are disjoint. 
From Lemma 13.11 we know that there exists T c > so that for T < T e the support of fii is 
contained in \xi{t) — e, Xi(t) + e\. 

Take i = 1, . . . ,p + q — 1, and fix the other measures fij, j ^ i. From (|2.12p . (|2.15|) we see 
that the measure Hi is the equilibrium measure in an effective external field 

/ X ^T^—\ d ^M + \Vi{x), (3.12) 

and it is also the minimizer if we restrict to measures with total mass tkU)M) that are supported 
on [xi(t) — e,Xi(t) + e]. On this interval the external field (|3.12p is strictly convex (we use that 
the measures \ij with j ^ i are supported outside [xi(t) — e, Xi(t) + e]). The convexity implies 
that the support of fa is an interval, see e.g. [301 Theorem IV 1.10]. 

(b) The real analyticity of the density of fa in the interior of the support follows from [T2] , 
since the effective external field (|3.12p is real analytic on the support of fa . 

The convexity of (|3.12p implies that the density of [ii does not vanish in the interior of the 
support, and has square root behavior at the endpoints, see e.g. [H Lemma 3.5]. □ 



3.2 Varying n 

In the situation of Theorem 12.41 we have for each finite n, the number rikj of paths going from 
cifc to 6;. The finite n transition numbers are 



and in the limit we have 



An) _ nk,i 

lim ti, = tki for k = 1, . . . ,p, I = 1, . . . , q. (3.13) 

n— >oc ' 

The equilibrium problem depends on the transition numbers by means of the normalizations 
(|2.17p . For a finite n, we use the equilibrium problem for a vector of measures (fix, . . . , /Up+g-i) 
that have total masses 

j dni =4"i),i(i)> for i = l,...,p + g- 1 (3.14) 

instead of ([2~T7)l . Then the minimizer (/Xj \ . . . , p^_ q _-y) will also depend on n. 

Because of Theorem 12.101 and (|3.13p . there is a To > 0, so that for every T < Tq, each is 
supported on an interval [a^ n ' ) , pf"] (depending on n) so that parts (a) and (b) of Theorem 1 2 . 1 Ol 
hold. As n — > oo, we have that /i[ n ' ) — > /i^ and 
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for every i. 

In the steepest descent analysis that follows we will fix a large enough n. We will work with 
the n-dependent measures fx^ and intervals [a[ n \ Pi ], but for ease of notation we will usually 
not write the superscript (n). We trust that this will not lead to confusion. However, we will 

write il n , . A property that will be used a number of times is that 

n ^ki = n M IS an integer. (3.15) 

4 Riemann surface, ^-functions, A-functions 
4.1 Variational conditions 

As said before, we take a large n and consider the vector equilibrium problem with normalization 
(|3.14p and we assume that T is sufficiently small so that Theorem 12. 101 applies. So the measure 
fa is supported on the interval [at, Pi]. 

The variational conditions associated with the vector equilibrium problem (12.15[) are as fol- 
lows. For each i — 1, . . . ,p + q — 1, there is a constant Li £ K so that 



— Li, x £ [ofj , Pi], % 
>L t , x£R\[a i ,(3 i \. ^' L) 



We use Fi to denote the Cauchy transform of the measure \Xi, 

/■ft i 

Fi{z) := / ^—dia(x), (4.2) 

for i = 1, . . . ,p + q — 1. The function is analytic on C \ [on, Pi]. By taking the derivative 

of (|4.1|) and using the Sokhotski-Plemelj formula it follows that 

- F ii+ (x) - F it _{x) -2j2ai,jFj(x) + ^'(i) = 0, i£ [cn.A]. (4-3) 

Lemma 4.1. TTie variational inequality (|4.ip is strict for x £ [f3i+i,on) U (/3,-,Oj_i], where we 
put ao = +oo emd /3 p + 9 = — oo. 

Proof. On both gaps (/3j+i, «i) and o^-i), the left-hand side of (|4.1[) is a real analytic function 
of a; whose first derivative is 

-2j2di,jF j (x) + ^Vf(x) (4.4) 



T 

3 



and whose second derivative is 



2 E a « ^ (^]2 + ^"(*)- ( 4 - 5 ) 

Each term in (|4.5[) is positive and so the left-hand side of (|4.1[) is strictly convex on both (Pi+i, ai) 
and {Pi, cti—i), which proves the lemma. □ 
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Figure 9: Four-sheeted Ricmann surface for the example in Figured] 



4.2 Riemann surface 

We construct a Riemann surface TZ as follows, compare with [17]. The Riemann surface has p + q 
sheets which we denote by IZj, j = 1, . . . ,p + q. Each sheet is associated with a vertex of the 
graph, that is, with either a starting point or an ending point 6;. We choose the numbering 
so that IZk is associated with ak for k = 1, . . . ,p, and 1Z p +i is associated with bi for I — 1, ... ,q. 

Recall that we use i as a label for the edges of the graph, and we write k = k(i) and I = 
if i labels the edge (ak, h). Then the p + q sheets arc defined as 

TZ k =C\ |J [oti,Pi], k = l,...,p, (4.6) 

i: k(i) — k 

TZ p+l =C\ |J l=l,...,q. (4.7) 

i:l(i)=l 

The sheets are connected as follows. For each i = 1, . . . ,p + q — 1 we have that Hk(i) is connected 
to sheet IZp+i^ along the interval [cii,Pi] in the usual crosswise manner. See Figure [5] for a 
picture of the Riemann surface for the example in Figure |H Note that the cuts [a^, Pi] are always 
between a sheet IZk which is among the first p sheets and a sheet 1Z p +i which is among the last 
q sheets. The cuts strictly move to the left if we go from one sheet to the next among the first 
p sheets, and similarly, among the last q sheets. 

The Riemann surface depends on n, since the endpoints and pf 1 depend on n. The £- 
and A-functions that we define from the Riemann surface will also depend on n. Wc will not 
indicate the n-dependence, as already mentioned before. 

The Riemann surface TZ has p + q sheets and 2(p + q — 1) simple branch points. Therefore 
by Hurwitz's formula (see e.g. [55]) its genus g satisfies 

2g-2 = -2(p + q) + 2(p + q-l) = -2 

so that g = 0. The fact that the genus is zero will be helpful in the construction of the global 
parametrix in Section [7.21 
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4.3 ^-functions 

We define the ^-functions as follows. Recall that Fi is given by 

&(«) = - F i( z ) + ^( z ~ a k), k = l,...,p, (4.8) 

i: k(i) — k 

tM*)= E Fi &- T (l-t) ( Z ~ bl) ' 1 = 1,...,?. (4.9) 

i:l(i)=l ^ ' 

We consider £j as an analytic function on the jth sheet IZj of the Riemann surface with a 
pole at infinity. Moreover, these functions define a global meromorphic function on TZ as the 
following result shows. 

Theorem 4.2. Consider the function £j{z) on the jth sheet IZj, j — 1, . . . ,p + q. Then these 
functions are compatible along the cuts [at, fa] of the Riemann surface in the sense that 

Zkw, + W = Zv+m,-(x), ie[ ft] (410) 

for every i = 1, . . . ,p + q — 1. Hence the £j -functions can be extended to a global meromorphic 
function £ defined on the Riemann surface TZ. 

Proof. Fix i G {1, . . . ,p + q — 1}. On the interval [at, Pi] we have by definition (|4.8[) (|4.9|) that 

£k(t),+ (a:) - Z p+ i(i),-(x) 

= - E f jM x )~ E p j>-W + Tta _ t) ( x -^- *) a *(i) - tb m)) 

j-.k(j)=k(i) r-i(j)=m [ ' 



-F it+ (x) - F i( _(x) - ( £ ^(s) I + ^V'{x), (4.11) 



V?¥* fc(i)=fe(i) or *0O=J(i) 

where we used the definition 13f> (|2. 14)) of Vi. 

If j 7^ i is such that = fc(i) or l(j) = l(i) then ajj = 1/2. For all other j ^ i we have 
dij = 0. Therefore 

\j¥=i- k(j)=k(i) or l(j)=l(i) J jjti 

Using this in (|4.11[) and recalling the variational equality f|4. 3[> . we see that that Cfc(i),+ ( a; ) — 
£p+l(i) ,-(x) for a; e [a l ,f3 l ]. 

The other equality £fc(i),-(a;) = £p+l(i),+ ( x ) follows in exactly the same way. □ 

In the following two lemmas we collect some more properties of the ^-functions that will be 
needed in what follows. 

Lemma 4.3. The ^-functions have the following behavior as z — > oo 
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Proof. From P~2"j) and (pHI)) it follows that 

(n) 

F, W _ifS + o(£)-^ + o(£) (4.14) 
asz^ cxi. Recall that we work with n-dcpcndcnt transition numbers. Since 

E,(n) _ ™fc \ ' (n) _ ^£ 

h(i),l(i) ~ n ' ^(O.lW - n 

i: k(i) — k i: 

the asymptotic behaviors (|4.12p and (|4.13p follow immediately from the definitions (|4.8p - (|4.9|) . 

□ 

Lemma 4.4. For any i = 1, . . . ,p + q — 1 we have 

Z k{f) (z)dz = -2Trit ( £l )tm , (4.15) 
f P +i(i)0) & = 27ritg), Ki) , (4.16) 
= 0, i^{*W.P + K0}. ( 4 -!7) 



where Ci denotes a counterclockwise contour surrounding the interval [aj,/3j] and not enclosing 
any point of the other intervals [aj,f3j], j =/= i. 

Proof. From (f^TT^j) and the definition of C, it follows that 

F t (z)dz = 27rHg ))i(i) 
and 

/ Fj(z)dz = iij^i. 
Jd 

The lemma then follows from the definitions (|4.8|) - ()4.9p . □ 

Remark 4.5. Our definition of the ^-functions differs slightly from the one used in [5]. If £j 
denote the ^-functions in [5] then we have 

~ z 

In the present form the formulae are more symmetric. 

4.4 A-functions 

We define the A., -functions as 

X k (z) = c k + f £ k (s)ds, k = l,...,p, (4.18) 



K+i{z) = Cp+i + £,p+i(s)ds, 1 = 1, ...,q, (4.19) 
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where 

ik '■= min{i | k(i) = k}, and ii := min{i | l(i) = I}, 

and where the path of integration in the integrals (|4.18p and (|4.19j) lies in C \ (— oo,/3j fc ) and 
C\ (—00, Pi ), respectively. The functions Xk{z) and A p +/(z) are defined with a branch cut along 
the intervals (— oo,Pi k ] and (— oo,/?jJ, respectively. 
We choose the constants Cj in in the following way. 

Lemma 4.6. We can (and do) choose the constants Cj in (I4.18p - (I4.19|) in such a way that 

Re(Afc (i) , + (A)) = Re(A p+;w , + (/J t )), (4.20) 

for every i — 1, . . . ,p + q — 1. 

Proof. We use the fact that the graph G — (V, E, t) is a tree. We can iteratively 'undress' this 
tree as follows. We start with G\ = G. Next we choose a leaf vertex v and set 

G 2 =G 1 \ v, 

i.e., G2 is the tree obtained by removing the leaf and its corresponding edge from the tree G\. 
We iteratively repeat this procedure and obtain in this way a chain of nested trees 

G = Gt D G 2 D ••• D G m , (4.21) 

where each Gi is obtained from G^_i by removing one leaf. Obviously the last non-empty tree 
G\v\ in this chain consists of a single vertex vj. 

We freely choose the corresponding constant Cj . Next we use induction on k = \V\ — 1, . . . ,2,1, 
to fill in all the remaining constants Cj so that each time (|4.20j) is satisfied. This is possible since 
Gi-i \ Gi consists of a single vertex Vj, which is a leaf of Gr»_i, and hence we have exactly one 
condition (|4.20p to fix the integration constant Cj of this leaf. Thus we see that the conditions 
(|4.20p can indeed be imposed on the constants Cj in (|4.18p - (l4.19p . □ 

Properties of the A-functions that we will need are stated in the following lemmas. The first 
is a reformulation of the variational conditions (|4.ip . 

Lemma 4.7. We have for i = 1, . . . ,p + q — 1, 

Re(A fc(i) , ± (*) - W ± (*)) { = J 11 j*'g' A] _ (4.22) 

Strict inequality in (|4.22p holds for x £ [ft+i, cti) U ((3i, cti-i] (where olq = +00 and (3p+ q = —00 ). 
Proof. Observe that by P~5]) - (|431) and (f4~THjl - (j4~T^]> we have that 

AfeW= Yl hog-^d^{x) + ^-{z-a k f + ~c k , (4.23) 



i: k(i)—k 



/l 1 
log— — dfa(x)- 2T(1 _ t) {z-bi) 2 + c p+l , (4.24) 

for certain real constants Cj, j — 1, . . . ,p + q. Therefore, for i = 1, . . . ,p + q — 1, 

Re (Afc (i) (z) - X p+ i(i)(z)) 

= 2 /log 1 dm (x)+ V / log 1 dfij (x) 

I \z — x\ ' / \z — x\ 

+ -Re^(z) + — (a fe(l) - 6 ;(l) ) 2 + c k(i) ~ c p+l[i) , (4.25) 
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where we used the definition (|2.14|) of Vi. Then by the variational conditions (|4.1[) we have 

Re(A feWi± (x) - \,+i(i),±{x)) > Li + —{a k{:i) - b l{t) ) 2 + c fc ( i) - c p+l(i) (4.26) 

for x E R with equality for x € [ccj, /?,•]. The constant in the right-hand side of (|4.26[) is equal to 
zero because the constants Cj are chosen so that (I4.20[) holds. 

The strict inequality for x € [A+i, on) U (fli, c^-i] is a consequence of Lemma l4.il □ 

Lemma 4.8. As z — * oo we have that 



M z ) = — (z - a k y logz + 4 + O - , fc = l,...,p, (4.27) 

2,1 1 n \ z J 

1 / 1 2 , mi , _ ~ ^ Z' 1 



A ^' (z) = ~ 2r(l-t) ( z — fy) + ~ log 2 + £p+i + Q . l = l,...,q. (4.28) 

Proof. This follows from the definitions 1|4.18[) and (|4.19p and the asymptotic behavior (|4.12p 
and P~T3f of the ^-functions. □ 

The next lemma will be a consequence of Lemma 14.41 

Lemma 4.9. For k = 1, . . . ,p, we have that 

exp(n(Ajk, + (a;)-Afc,-(x))) = l, /or i6l\ |J [a,, ft]. (4.29) 

z: k{i) — k 

For I = 1, . . . , q, we have that 

exp(n(X p+ i t+ (x) - Xp+i-(x))) = 1, /orireR\ (J [a^]. (4.30) 

i: l(i)=J 

Proof. Fix fc = 1, . . . ,p and let a; G R \ (J^. fe ^^_ fc [aj, /3j]. By definition of Xk and by contour 
deformation we have that 

Xk,+{x) - Xk-(x) = J £k(z)dz 

where C is some closed contour surrounding some of the intervals [aj,/3j]. From Lemma 14.41 
it follows that each of the enclosed intervals gives a contribution to the integral of the form 
±27ritj|."\ Since each t^j is a rational number with denominator n, indeed t^f} — rikj/n, we 
conclude that n(Xk,+ (x) — (x)) is a multiple of 2-7ri, and so we obtain (|4.29j) . 

The proof of ([430]) is similar. □ 

Lemma 14.91 will be important for the steepest descent analysis in Section O 



5 First transformation of the RH problem: Normalization 
at infinity 

In the following sections we describe the steepest descent analysis FnX^T^^i-iJiof the 
RH problem 12. Ill Throughout the steepest descent analysis the following simple lemma will be 
repeatedly used. 
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Lemma 5.1. Assume that the matrix function Y{z) satisfies the jump condition Y + (x) = 
Y^(x)J(x) for Let A(z) and B(z) be matrix functions with A{z) entire and B{z) analytic 

mC\R. Then 

X(z) := A(z)Y(z)B(z) 

satisfies the jump condition 

X+{x) = X_(x) (BZ 1 (x)J(x)B + (x)) . 

The point will be to choose appropriate transformation matrices A(z) and B{z) in order to 
bring the RH problem 12.111 to a simple form. 

Let Y be the solution to the original RH problem 12.111 The first transformation Y i— ► X 
serves to normalize the RH problem at infinity. To this end we define X — X(z) as 

X(z) = L- n Y(z)D(z) n , (5.1) 

where we define the diagonal matrices 

1 



D(z) = diag ( ( exp ( X k (z) - -^(z - a k f 



k=l 

1 



^{^W + wTQ^ JjJ' (5 ' 2) 

and 

L = diag(exp(ci), . . . , exp(c p+q )), (5.3) 

where the constants cj are as in (I4.27j> — ()4.28|) . From Lemma f5.1l it follows by straightforward 
calculations that X = X(z) satisfies the following RH problem. 

RH problem 5.2. 

(1) X is analytic in C \ K; 

(2) For x € K, we have that 



where the blocks Ji t t, Ji,2 and J2.2 (of sizes p x p, p x q and q x q, respectively) have the 
following form. 

(a) >A,2 is a full p x q matrix with entries 

{Ji,2{x)) kl = exp(n(X p+ i t+ (x) - X k ,-{x))) , xeR. (5.5) 

(b) Outside of the intervals [ai,/3i], Ji ; i and Ji,2 are identity matrices 

p+q-l 

J 1 , 1 (x)=I p , J 2 ,2{x)=I q , a;6l\ (J [cn,pi\. (5.6) 

8=1 
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(c) On the interval [aj, /?,•], J\ t \ and J 2j2 are diagonal matrices with ones on the diagonal, 
except for 

( J i,i( a; ))fe( l ),fc( l ) = ex P i n (h(i),+ { x ) ~ -( x ))) , x ^ (atuPi), (5-7) 



(3) As z — > oo, we have that 



X(z)=I p+q +0(l/z). 



(5.9) 



Here we used Lemma 14.91 to see that the diagonal entries of (|5.6|) , as well as most of the 
diagonal entries of (|5.7p - (|5.8p are equal to 1. On the other hand, we used the asymptotic 
conditions (|4.27j) - l|4.28p and (|5.ip - (|5.3[) to see that the RH problem for X is normalized at 
infinity in the sense of (|5.9[) . 

Let us illustrate the jump matrices for the example with p = q = 2 as in Figures 0] and [HI In 
that case the jump conditions are written as 



X+ = X. 



on the interval [ai,/?i], 



X + = X- 
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exp(n(A4 ! + — 
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on the interval [a 3 ,/3 3 ], and 



X+ = X- 



(l exp(n(A 3 - Ai)) exp(n(A 4 - Ai))^ 

1 exp(n(A 3 - A 2 )) exp(n(A 4 - A 2 )) 

1 

\0 1 / 



onR\uLiK&: 



6 Second transformation of the RH problem: Gaussian 
elimination and opening of global lenses 

— Afc(i) _ is a purely imaginary constant. So the 



On the interval [cti,f3i] we have that Ap+j^) . 
entry 



(6.1) 
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is constant with modulus one. It would be an ideal situation if, except for (|6.ip , all entries of 
the matrix Ji^ix) for x S R, are exponentially decaying as n — * oo. That would happen if for 
every k = 1, . . . , p, and £ = 1, . . . , q, we have 



(a) if 4™ } = then 



(b) if 4" } > then 



where i — k + I — 1 . 



Re A p+ ; i+ (x) < ReAfc i _(a;), 
ReA p+ ; i+ (a;) < ReA fe: _(a;), x eR \ [a i} (3i], 



However, this will not happen in general, and that is why we need a second transformation 
1 h T in the steepest descent analysis of the RH problem in which a number of unwanted 
entries of the jump matrices arc eliminated. In particular those entries of J\^{x) that could 
potentially be exponentially increasing as n — > oo. 

To this end we will open global lenses [3] and apply Gaussian elimination inside each of the 
lenses. This construction will be systematic and may have interest in its own right. The proof 
that appropriate lenses exist will be a consequence of the maximum principle for subharmonic 
functions. 

The opening of global lenses can be conveniently described in terms of the right-down path 
in Proposition ^. 21 We start at the top left entry (1, 1) of the matrix of transition numbers (t^j) 
and walk along the path until we arrive at the bottom right entry (p,q). During this walk, we 
will open global lenses in an appropriate way. The precise action to perform depends on whether 
we are taking a vertical (down) or a horizontal (right) step along the path. 

6.1 Construction of global lenses: Vertical step 

First we construct the global lenses for a vertical step along the lattice path. Thus assume that 
i € {1, . . . ,p + q — 2} is such that 

(k(i),i(i)) = (k,i), (6.2) 

(*(i + l),l(t + l)) = (A + 1,0, (6-3) 

for certain k = 1, . . . ,p — 1 and I = 1, . . . , q. 

From (|4.27p we obtain the asymptotic behavior 

Xk+i(z) - A fe (z) =^( fl fe - a k +i)z + 0(\og \z\), (6.4) 

as z — * oo. 

We also note that Re A^ and Re A^+i are well-defined and continuous on C. Indeed, we have 
by that 

r .KU)=K V " xl 

log rTj ^j(^) + 7ffi Re ( z - afc+i) 2 + c fc +i, (6.6) 



for certain constants c k and Ck+i- 

The representations (|6.5[) and (|6.6p also show the following. 
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Lemma 6.1. The function z i— » Re (Afe+i (z) — Afc(z)) is superharmonic on C \ (J ., ^.ij [ay, 
and subharmonic on C \ (J^. fc (j) =fc+1 [atj , /3j] • 

Proof. It is a standard fact from potential theory that any function of the form z i— > J log uz^i d/j,(x), 
with /i a positive measure with compact support, is superharmonic on C and harmonic on 
C\supp(^i), see e.g. [SOI Chapter 0]. Thus by (|6.5p and (|6.6|) we have that z i— > Re Afc+i(z) is su- 
perharmonic on C and harmonic on C \ (Jj. fc(j)=fc+i [ a i > /%]> while z i— > —Re (z) is subharmonic 
on C and harmonic on C \ (J^. k{j)=k \ a jifij\ 

Since the two sets Uj- fc(j)=fc[ Q! i) /%] an d Uj- fc(j)=fc+i ft] are disjoint, the lemma follows. 

□ 

We next define the open sets 17+, 17 _ C C as follows 

17+ := {z g C | Re (A fc+1 (z) - A fe (z)) > 0} (6.7) 

17_ := {z G C | Re (A fc+1 (z) - A fe (z)) < 0}. (6.8) 

We also denote 

17+, oo := {z G 17+ | 3 a connected path in 17+ from z to oo} (6-9) 

17_ j00 := {z G 17_ | 3 a connected path in 17_ from z to oo}. (6.10) 

In other words, 17+ !00 is the union of the unbounded connected components of 17+, and similarly 
for S7_ !00 . 

The open sets 17+ !00 , 17_. 00 satisfy the following properties. 
Lemma 6.2. (a) For each e > there exists R > so that 

{z G C | |z| > R, — k/2 + e < argz < tt/2 - e} C 17+ i00 , (6.11) 

and 

{z G C | |z| > R, ir/2 + e < argz < 3tt/2 - e} C 17_ i00 . (6.12) 

In particular, 17+ iOC Zzes £o the right o/17_ 00 . 

(b) Bo£/i 17+ !OC and 17_ j00 are connected. 

Proof. Part (a) follows from (|6.4p and the the fact that a k > afc+i- Part (b) follows from part 
(a) in a similar way as in Proof of Lemma 2.4], to which we refer for further details. □ 

Lemma 6.3. We have 

ol% G 17+ :00 and i+1 G S7_ !00 , (6.13) 
where we recall i is related to k as in (|6.2[) and (|6.3p . 

Proof. By applying the variational conditions (|4.22p twice, first with the index i and then with 
i + 1, we obtain 

Re(A fc +i(x) - X k (x)) = Re(A fe +i(a;) - \ p +i(x)) > 0, at e (6.14) 

and the inequality (|6.14p is strict for a; = on because of the statement about the strict inequality 
in Lemma 14.71 Hence a, G 17+, and in a similar way we obtain /3,+i G 17_. 
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To show that on belongs to the unbounded component of f2 + we argue as in (TBI Proof 
of Lemma 2.4]. What we use is that \J j: k ( j)=k [ay, 0j] lies to the right of [jj. k ( j)=k+1 [a 3 ■, f)j], 
and that a.i is the left-most point of (J .. k (j)=k [ a j> @i\ > an( ^ that A+i ^ s ^ nc right-most point of 

Uj:fc(j)=jfe+l[ a Ji A']- 

The proof is by contradiction. Suppose that ctj does not belong to the unbounded component 
of f2+. Then the set 

f2 +jQi := {z G fl + | 3 a connected path in 17 + from z to a{\ (6.15) 

is bounded, it is symmetric with respect to the real line, and it contains a,. Also Re(Afc + i — 
Afe) is zero on the boundary of Q+, ai and strictly positive inside of by construction. 

Since subharmonic functions satisfy a maximum principle, it follows that Re (Afc+i — X k ) is not 
subharmonic on all of f2+ jQi . Then by Lemma 16.11 we conclude that Q+, ai has a nonempty 
intersection with Uj- fc(j)=fc+i [ a i> ft] ■ ^ny point of intersection lies strictly to the left of Pi+i, 
since /3i + i G f2_ and f3i+\ is the right-most point of (J .. fc ^-)_ fc+1 [ctj, /%]- Then, because of 
symmetry in the real axis, it follows that £l+. Qi surrounds the point ft+i. The set 

r2_^ i+1 := {z G 0_ | 3 a connected path in fL_ from z to (6.16) 

is then bounded, and it does not intersect with \Jj. k ij\ =k .[oij, J3j], By Lemma 16. 1[ we then 
have that Re(Afc+i — Afc) is superharmonic on f2_ j g 1+1 . However, Re(Afc+i — X k ) is zero on the 
boundary, and strictly negative inside of Q— t p i+1 , which gives a contradiction with the minimum 
principle for superharmonic functions. 

Thus a.i belongs to the unbounded component of f2+, and similarly Pi+i belongs to the 
unbounded component of Q_ . □ 

Let Xq > be such that 

Re(A k /(x) - Xp+i'(x)) > 0, 

for all x G (— oo, — Xq) U (Xq, oo) and all k' = 1, . . . ,p, I' — 1, . . . , q. The existence of such a 
constant X follows from (|4~271l - (|4?28| . 

The next result is an immediate consequence of Lemma 16.31 

Theorem 6.4. There exist two simple closed contours 

T +ti C fi+,00 and T- :l C H-.oo 

such that 

(a) L +! i surrounds the interval 

(b) surrounds the interval [a p + q _i, /3j+i]. 

(c) ooi/i and L +i i intersect the interval m exactly one point, which we denote 
by Xi, yi, respectively. We have Xi < yi. 

(d) both L_^ and T+i Ziawe one extra intersection point with the real axis, which lies inside the 
interval (-co, —Xq), (Xq,oo), respectively. 

As in fSjj, the contours T+ t i and F-^ will be called global lenses. 

Remark 6.5. Instead of taking closed contours, one might also take and to be un- 
bounded, tending to infinity in the right and left half of the complex plane, respectively, and 
both intersecting the real line in exactly one point in the line segment (/3i+i, on). This is the 
construction that was used in [IB] , 

An illustration of Theorem 16.41 is shown in Figure [TUJ When convenient we also define 
xq = X , y = oo, x p+q -i = -oo and y p+q -i = -Xq. 
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Figure 10: The figure shows how to open global lenses between the two intervals [oij+i, 
(left) and (right). 



6.2 Construction of global lenses: Horizontal step 

Next we construct the global lenses for a horizontal step along the lattice path. Thus assume 
that i € {1, . . . ,p + q — 2} is such that 

(fc(i),J(*)) = (M). (6-17) 
(fc(i + l),J(i + l)) = (M + l), (6-18) 

for certain k — 1, . . . ,p and I = 1, . . . , q — 1. 

The analysis of Section 16.11 can be adapted to the present case. Let us discuss the main 
points. We now define the open sets f2+ and fi_ as follows 

n+:={zeC\ Re (A p+; (z) - X p +i+i(z)) > 0} (6.19) 
0_ := {z e C | Re (A p+i (z) - A p+;+1 (z)) < 0}. (6.20) 

The unbounded regions 0+,oo) ^-,°o c ^ are a g am defined as in (|6.9j) ~ (|6. 10|) . Then Lemma [6721 
remains valid. 

Using the above definitions, Lemma l6.3l and Theorem 16 .41 both remain valid as well. Thus we 
can construct the global lenses r+ ^ and r_.j in exactly the same way as before. 



6.3 Gaussian elimination step 

Now we will show how to apply Gaussian elimination inside the global lenses. First, it might be 
worthwhile to recall the basic idea of Gaussian elimination in our present context, see e.g. [2"U] . 
Let 




(6.21) 



be a rank-one matrix of size p by q. Assume that one multiplies J on the left with the matrix 

I P - — e fe+1 e£. (6.22) 
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Here and below we use to denote the column vector with all entries equal to zero, except for 
the kth entry which equals 1. The length of will be clear from the context. Note that the 
outer product ek+iej in ()6.22f> is the matrix with all zero entries except for the (k + l,k) entry 
which equals one. 

The multiplication with (|6.22[) on the left is equivalent to applying an elementary row oper- 
ation to the rows of J, where to row k + 1 one adds — times row k. This row operation is 
such that the entries of row k + 1 of J are eliminated. 

Similarly assume that one multiplies J on the right with the matrix 

/ ? - ^±i e; ef +1 . (6.23) 

This is equivalent to applying an elementary column operation to the columns of J, where to 
column I + 1 one adds — times column I. This column operation is such that the entries of 
column Z + 1 of J are eliminated. 

Let us see how we can apply these ideas in the present context. The role of the rank-one 
matrix (|6.21j) of size p by q will be played by the top right submatrix J\^(x) of the jump matrix 
(|5.4[) . cf. ()5.5|) . The mechanism to multiply the jump matrix on the left or on the right with a 
transformation matrix of the form (|6.22[) - (|6.23p is to define a new RH matrix T(z) = X(z)B(z) 
for a suitable transformation matrix B(z) and subsequently apply Lemma 15.11 

Now we are ready to describe the Gaussian elimination in detail. This will be the next 
transformation X i— > T in the steepest descent analysis of the RH problem. 

Algorithm 6.6. (The transformation X i— > T) 

1. (Initialization.) We initialize T(z) := X{z). 

2. (Forward sweep.) For each i = 1, . . . ,p + q — 2 we open the global lens T—j in Theorem \6.4\ 
and we update, in case of a vertical step (|6.2p - (|6.3p , 



, v _ f T(z) (l p+q + cxp(n(A fe+ i(z) - \ k (z)))e k el +1 ) , inside the lens T_ 
1 T(z), elsewhere, 



T 

and in case of a horizontal step (|6.17p - (|6.18p . 



T 



r z \ = | T ( z ) (lp+q - exp(n(A p+; (z) - X p+ i +1 (z)))e p+l+1 eJ + ^ 

I T(z), elsewhere 



inside 



(Backward sweep.) For each i = p + q — 2, ... ,2,1 we open the global lens in Theorem 
\6.4\ and we update, in case of a vertical step (16.2p - ()6.3|) . 

inside the lens 
elsewhere, 



T(,\ - / T(z ) ( 7 p+9 + exp(n(A fc (z) - X k+1 (z)))e k+1 ef) 
[ >~ { T(z), 

and in case of a horizontal step (|6.17p - (|6.18p . 

T(z) = 



T(z) (lp+q - exp(7i(A p+ /+i(z) - X p+ i(z)))e p+ ie^ +l+1 ^j , inside T + . t 
T(z), elsewhere. 



Incidently, we note that the forward and backward sweeps in the above algorithm commute. 
But one is not allowed to change the order in which the index i varies inside the sweeps. 
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It is easy to see that Algorithm 16.61 does not change the jump matrix in (|5.4p except for its 
top right submatrix Ji^ix). To see what happens with the latter, we have to distinguish between 
different regions of the complex plane. 

First assume that x belongs to one of the intervals (jjj, Xj—i) D [ay, j = 1, . . . ,p + q — 1. 
(Recall that y p + q -i = —Xq and xq = Xo.) From Theorem l6.4l we see that x lies inside the global 
lens r_^ precisely when i = 1,. .. ,j — 1. Hence the 'relevant' indices in the forward sweep in 
Algorithm 16.61 are i — 1, ... , j — 1. During the corresponding operations, the entries in rows 
1,2,..., k(j) — 1 and columns 1,2, ... , l(j) — 1 of the matrix Ji^(x) are cancelled by Gaussian 
elimination. 

On the other hand, Theorem 16.41 shows that x lies inside the global lens r + i precisely when 
i = p + q — 2,p + q—3,...,j. Hence the relevant indices in the backward sweep in Algorithm l6.6l 
are i = p + q — 2,p + q — 3, . . . , j. During the corresponding operations, the entries in rows 
p,p — 1, . . . , k(j) + 1 and columns q,q — 1, . . . , l(J) + 1 of J\.2{x) are cancelled by Gaussian 
elimination. 

It follows that at the end of the two sweeps in Algorithm 16. 6[ all the entries of the rank-one 
matrix Ji^ix) are eliminated, except for the (k(j),l(j)) entry which equals 



Recall that in the above description, we assumed that x £ (yj,Xj_i) D [aj,&j]. Next, let us 
assume that x belongs to one of the gaps {xj,y,j) for certain j. We can then repeat the above 
arguments and find that at the end of Algorithm 16. 6[ all the entries of Ji^(x) are eliminated 
except for two of them. In case of a vertical step (|6.2)) — (|6.3|) these are the (k(j),l(j)) and 
(k(j) + 1, entries, which are given by 



respectively. But by the variational inequality in (|4.22p , which is strict according to Lemma 14. 1[ 
we see that both entries are exponentially small for n — > oo. A similar argument applies in case 
of a horizontal step (|fTT7|) - ([6Tl| . 

Finally we should note that by the operations in Algorithm l6.61 the RH matrix T{z) also has 
a jump on each of the contours r +ji and For example, in case of a vertical step (|6.2p ~ (|6.3p 

the jump matrix on the contour takes the form 



see Algorithm 16.61 But by our assumption that ^ C f)+ we see that this jump matrix is 
uniformly exponentially close to the identity matrix when n — > oo. A similar argument holds for 
the jumps along the contours r_^ C f2-. 

Summarizing, we established that T satisfies the following RH problem. 

RH problem 6.7. 



exp(n(X p+l{jh+ (x) - A fe(i) -(x))). 



exp(n(X p+ i U} {x) - \ k (j)(x))), exp(n(\ p+lU} (x) - \ k{j)+1 (x))), 



I p+q ± exp(n(A fc (z) - A fe+ i (z)))e k+1 el, 



(1) T is analytic on C \ (M U (J 



IP+9-2 



(r+,,uL,,)). 



(2) For xeRUlJ: 



IP+9-2 



(r +i i U r_^)) we have that 



T + {x)=T_{x)J T (x) 



(6.24) 



where Jt{x) satisfies the following 
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Figure 11: The figure shows the contours in the RH problem for the matrix T for the example of 
Figures 2] and [5] We have three intervals [0*,/%], i — 1,2,3, and four global lenses T±,i, i= 1,2. 



(a) For x £ (j/j,a;,_i) D [aj,/3j], i — 1, . . . ,p+q—l, we have that Jt{x) equals the identity 
matrix, except for the 2x2 block lying on the intersection of rows and columns k(i) 
and p + l(i), which is given by 



exp(n(A fc(4 ) !+ (a;) - A fc(i) _(x))) exp(n(X p+ i(i) t+ (x) - A fe (;) _(x))) 
exp(n(\ p+l{lh+ (x) - Ap + j W) _(x))) 



(6.25) 



(b) For x G ( — oo, Vp+q—i) U (UiO^iff*)) U (xo,oo) we have that Jt{x) is exponentially 
close to the identity matrix as n — > oo. 6o</i uniformly as well as in L? sense. 

(c) For x G (J i (r+ i j U r^^), the jump matrix Jt is also exponentially close to the identity 
matrix as n — > oo m uniform sense (and therefore also in L? since the contours T± t i 
are compact). 



(3) As z —f oo, we have that 



T(z) = I p+q + 0(l/z). 



(6.26) 



Let us illustrate this RH problem for the example in Figures [4] and [9] Then we have three 
intervals [aj,/3j], i — 1,2,3, and four global lenses between them: see Figure [TT1 
The jump conditions (|6.24p - (|6.25p can now be written as 



/exp(n(Ai. + - Ai _)) exp(n(A 3 ,+ - A x _)) 0\ 



V 



1 o 

exp(n(A 3 .+ 






A 3 ,-)) 
V 



on the interval (yi,Xo) 3 [c*i)/?i] 


















exp(n(A 2 .+ - 


- A 2 ,_)) exp(n(A 3:+ 


-A a ,_)) 











exp(n(A 3 .+ 


-A 3 ,-)) 















V 



32 









Vi 


ok 


/fa k-i K 



Figure 12: A local lens around the interval [aj,ft]. 
on the interval (j/2,^1) 3 [a2,ft>], and 



T+ = T_ 



/l ^ 

exp(n(A 2 ,+ - A 2 ,_)) exp(n(A4. + - A 2 ._)) 

10 

\0 exp(n(A 4 ,+ - A 4 ,_))7 



on the interval (j/3, £ 2 ) D [03, /3a]. The jump matrices on the remaining contours (—00,2/3), 
(^2,2/2)5 (xi,Ui), (xo,oo), r +i i, r_ i i, r + 2 and r_ j2 in Figure [TT] are all exponentially close to 
the identity matrix as n —* 00. 

7 Final transformations of the RH problem 

The jump matrices Jt in the RH problem for T are nontrivial only in the 2x2 block given by 
(|6.25|) on the interval [o!j,ft]. 



7.1 Third transformation: Opening of the local lenses 

In the transformation T 1— ► S of the RH problem we transform the oscillatory entries of the jump 
matrix (|6.24[) - (|6.25p along each interval [c*i, ft] into exponentially decaying ones. To this end we 
open a local lens around the interval [o^, ft] = [on, ft]. Since the RH problem is locally of size 2 
by 2 this can be done in the standard way [HJ [TS] . 

For each i — 1, . . . ,p + q — 1 we open a lens around the interval [at, ft], i = 1, . . . ,p + q — 1 
as in Figure [12] so that 

Re (A fe (j) - Xp+i(i)) < 

on the lips of the lens. It follows from an argument based on the Cauchy-Riemann equations that 
this is indeed possible, cf. [TT]. The lenses are small (so called local lenses), and in particular 
they do not intersect with the global lenses T±^. We use and to denote the upper and 
lower lip of the lens, respectively. 
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Wo define the matrix function S as follows 



T(z) (l p+q - exp(n(A fc(l) (z) - \ p+l{i) {z)))e p+m el (i) 

in the upper part of the lens around [aj,/3j], 

S ( z ) = < T{z) (l p+q + exp(n(A feW (z) - X p+l(t) {z)))e p+l(i) el {i) 

in the lower part of the lens around [on, Pi] 
T(z), outside of all the lenses. 

Then S satisfies the following RH problem. 
RH problem 7.1. 

(1) S is analytic mC\(RU U 4 ( r +,* U T- ti ) U U i ( i +,i u L -,i))'> 

(2) For igKU Ui( r +,i u r -,0 u UiC-^+.t u we have that 

Sj r {x) = S-(x)Js(x) 

where Js(x) satisfies the following 



(7.1) 



(7.2) 



(a) For x £ [yi, Xi—x] U U L-.i we /lave i/iai Js(x) is the identity matrix except for the 
2x2 block on the intersection of rows and columns k(i) and p + l{i), which is given 
by 

' o r, ,. , 



-1 

1 (n 

exp(n(A fc(l) - lj ' 



/or x e L±,i, 



and &?/ 



exp(n(A fe (j) i+ (x) - A fc(i) _(x))) exp(n(A p+iW!+ (:z;) - X k ^_(x))) 
exp(n(Ap +iW , + (x) - A p+iW _(x))) 



(7.3) 
(7.4) 

(7.5) 



/or x € [yi,aj] U [/3 4 , x 4 _i] . 

(b) On </ie other parts of the contour we have Js(%) — Jt{x) and Js(x) is exponentially 
close to the identity matrix as n — > oo, both uniformly and in the L 2 sense. 



(3) As z — > oo, we /iave i/iai 



5(z) = I, 



P+9 



•0(1/*) 



(7.6) 



From standard arguments based on the Cauchy-Riemann conditions |llj it follows that the 
local lenses Li can be chosen so that the jumps on Li in (|7.4p are uniformly exponentially close 
to the identity matrix, away for a neighborhood of the endpoints on, fy. 



7.2 Global parametrix 

In this subsection we build a global parametrix p(°°\z), which will be a good approximation to 
the RH problem away from the endpoints on, Pi, i — 1, . . . ,p + q — 1. The construction will be 
quite similar to the one in [lOj . 

We will construct the matrix function pt°°) (z) such that it satisfies the following RH problem, 
obtained from RH problem 17. II by ignoring all exponentially small entries of the jump matrices. 
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RH problem 7.2. 

(1) P(°°)(z) is analytic in C\ Ufi 9-1 ; 

(2) For x S Ui( a i' @i)> we have that 



p| co) (x)=P H (x)J PW (x) 



(7.7) 



where the jump matrix Jp(oo> (.t) equals 



Jp(oo) (x) 

/or a; € (aj,/3i). 
(3) As z — > oo, we /iawe i/iai 



/ h(i)-l 





V o 







I p +i(i)-i-k(i) 

-1 





o \ 






(7.8) 



Iq-l({)J 



P(°°\z)=l p+q + 0{l/z). 



To solve this RH problem, we will use the fact that the Ricmann surface 1Z in Section | 
has genus zero. From general algebraic geometry [26j . this implies the existence of a rational 
parametrization 



£ = z = z{v) 



(7.9) 



where v runs through the extended complex plane C (Riemann sphere). 

The w-plane is then partitioned into p + q disjoint open sets Qj, j — 1, . . . ,p + q, where 
Slj is defined as the inverse image under (|7.9[) of the jth sheet TZj of the Riemann surface. 
Correspondingly we have p + q inverse functions Vj(z) of (|7.9p such that 



j = l,...,p + q, 



(7.10) 



is a bijection. We use Vj(po) to denote the image under this map of the point at infinity of the 
jth sheet TZj, j = 1, . . . ,p + q. Hence Vj(oo) <E flj. 

For i = 1, . . . ,p + q — 1, the common boundary of fifc(i) and fl p+ i^ in the w-plane, is an 
analytic curve Ci with a natural partition 



(7.11) 



where C+.i is the image of the upper side of the cut [ai,/3i] under the mapping v^, and is 
the image of the lower side. The two parts C± i meet at two points 7^ and 7. , which are the 
images of the endpoints cti and Pi of the cut [aj,/3i], respectively. 
Define the polynomial 



p+q-l 



g(v)= n ((«-7, (1) )(«-7f 2) )) 



(7.12) 



and its square root 
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which is defined as an analytic function in the u-plane, with a cut along the disjoint union of 
arcs U i C+,j. We assume yf g(v) ~ as v — > oo. 

We then construct a global parametrix P(°°\z) as in [10] . We define for z G C \ UilP^i A]i 

p(-)(.) = , /*(«) = ^ 7 - 13 ) 

where is the Lagrange interpolation polynomial for the points V\(oo), . . . v p + q (oo). That is, U 
is a polynomial of degree p + q — I so that 

/»(«i(oo)) = <5 ij7 j = l,...,p + g. 

The fact that P(°°'(z) in (]7. 13j) satisfies conditions (1) and (3) in RH problem l7.2l is immediate. 
For the jump condition (2) we need to show that 

/,(«*«,+ (»)) = ~fi(v p+lii) (*)) 1 

/iK+iw.+W) = /iKi),-( ji )). / 



These relations reduce to 



A+(«) = ^ v e C+,i, (7.15) 

fi,+ (v) = fi,-(v), torveC^i (7.16) 



and these jumps follow from (|7.13|) . since we have chosen the square root in yj g(v) with a cut 
along the union of arcs Ui^-+,i- 

7.3 Local parametrices around the endpoints: Airy parametrices 

In a small disk around the endpoints on and Pi of the interval [a!j,/3j] we construct a local 
parametrix p( Alr y) (2) involving Airy functions. Since the RH problem is locally of size 2x2 and 
the equilibrium measures all vanish as a square root, this can be done in the standard way [llj . 
We omit the details. 

7.4 Final transformation of the RH problem 

Using the global parametrix p(°°) of Section 17751 and the local parametrices p( Alr y) of Section 17731 
we define the final transformation S 1— ► R of the RH problem by 

T3( \ _ f S , (z)(P (Airy:) ) _1 (z), in the disks around on,j3i, i = 1, ...,p+ q - 1, C7 171 
[Z> ~ { 5(z)(P(°°))- 1 (z), elsewhere. 1 ' 

From the construction of the parametrices it then follows that R satisfies the following RH 
problem. 

RH problem 7.3. 

(1) R(z) is analytic in<C \ Y,r where is the contour shown in Figure [7731 

(2) R has jumps P+ = R-Jr on Xr, where 

Jr(z) = Ip+q + 0(1 /n), on the boundaries of the disks, 
Jr{z) = Ip+ q + <3(e~ cn(W+1) ), on the other parts of£ R , 
for some constant c > 0. 



3G 



r-,i r+, 2 




Figure 13: The figure shows the contours in the RH problem for the final matrix R(z) for the 
example of Figures 0] and [5J 



(3) R{z) = Ip+q + 0(1/2) as z — > 00. 

As n — > 00, the jump matrix Jr tends to the identity matrix both in L°°(Sfl) and in L (£r). 
Then as in [TTJ [T31 [T3] we may conclude that 

*<*>=^ +o G^+i)) (? - i8) 

as n -» 00, uniformly for z in the complex plane. This completes the RH steepest descent 
analysis. 



8 Proof of Theorem 12.4 



Now we are ready to prove the main Theorem 12.41 by unfolding the transformations of the RH 
steepest descent analysis. Compare with the proofs in the earlier papers [51 151 HP]. 
For a finite n we define the function p( n > as 



P 



W(x) = ~lm^ + (x), xe[a[ n \^\ i = l,...,p+q-l, (8.1) 



p+q-l 



p<- n \x)=0, xeR\ |J [a\ n) ,pl n) ]. (8.2) 



Here we write + (x) and a\ n \ /?[ to emphasize the rt-dependence. 



We recall that 



Im£ (n) --Im£ (n) --Im£ (n) = Imf (n) 



on the interval [a|™ , /?j ], so one has in fact several equivalent ways of expressing (|8.1 
By (|4.8|) . I|8.1[l and the Stieltjes-Perron inversion formula one has that 
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As n — > oo, we have that a\ — > on, /?■ — > fa and 



where 



lim p ( - n '{x) = pi(x), x € (at, fa), 

n — >oo 

dfj,i(x) 



(8.3) 



dx 



x e [ai,fa] 



is the density of the ith component p,i of the minimizer (/ii, . . . , p p+q -i) of the vector equilibrium 
problem with transition numbers {fk,l)- 

Now we show that the pi give indeed the limiting distribution of the non-intersecting Brownian 
motions. To this end we will use (|2.37|) . We start with the expression for the correlation kernel 
(|2.36p . which we restate here for convenience 



K n (x,y) 



2tti(x — y) 



(0 ••• W2,i(y) ••• w 2 , q (y))Y+\y)Y + (x) 



wi, p (x) 




From the first transformation Y i— > X in (15. ip 
n-dependence in the A-functions) 



V o J 

we get (we do not explicitly write the 



K n (x,y) 



1 



2ni(x — y) 



(0 ••• e n Vn.+to) ... e »Vh,,+(2/)) X^ 1 (y)X + (x) 



/ e -nXi 1+ (x)\ 




V o J 

From the second transformation I^Twe obtain for x,y S {otf 1 ^ , ft\ ) , 
From the third transformation T i— > S* in (|7.ip we get 



f T ,,1 — - ( - e nX k(i),+(y) e T, , + P riA p+i( l ).+(^) e T „ \ 

A " 11 ' yJ " 27ri(a; -y)\ k(l) p+l(t) 



x S?(y)S+(x) (e- nX ^+Me k(i) + e-" A *+«^+^e p+ *w) , (8.4) 



for i,i/6 ("^i i A' ) • Defining the function h n by 



(n) o{nh 



h n (x) := -Re(A fc(i ) !+ (x)) = -Re (A p+ j Wj+ (ie)), a; € [a,- ,# 
we see that (|8.4[) can be rewritten as 

„n.(/i ra (x)— h„(y)) 



(8.5) 



K n {x,y) 



2i:i{x — y) 



x S?(y)S+(x) (e-™ Im ( A M 1 ), + (-)) efeW +e" iIm ( A M I ),+ (-))e p+;(l) ) , (8.6) 
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tor x,ye {a™,$ n >). 

Now from (|7.18p it follows by standard arguments (e.g. Section 9]) that 

S^ 1 (y)S(x) = I + 0(x - y), as y -> x 

uniformly in n. Hence (|8.6|) takes the form 

K n (x, y) = e«(M*)-M,)) f S in(nIm(A fcW + (,)-A fc(l , + (,))) + A 

for x, y G (cti,f3i), where the O(l) term holds uniformly in n. Then by letting y — > a; and using 
l'Hopital's rule we find 

^„(ar,x) = ^(lm(cg )>+ ( a; )))+0(l) ) 
for x € (aj,/3j), or equivalently 

# n (z,a;)=np"(z)+0(i), 
by virtue of (jlTTj) . It follows from that 

lim -K n (x,x) = pi(x), xe(ai,(3i). (8.8) 

In a similar way one can prove that 

lim -KJx,x) =0, scGR\l Jfaj./Sjl. (8.9) 

This completes the proof of Theorem 12.41 
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